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Abstract 

In this paper we study direct and inverse problems for discrete 
and continuous time skew-selfadjoint Dirac systems with rectangular 
(possibly non-square) pseudo-exponential potentials. For such a sys¬ 
tem the Weyl function is a strictly proper rational rectangular matrix 
function and any strictly proper rational matrix function appears in 
this way. In fact, extending earlier results, given a strictly proper ra¬ 
tional matrix function we present an explicit procedure to recover the 
corresponding potential using techniques from mathematical system 
and control theory. We also introduce and study a nonlinear gen¬ 
eralized discrete Heisenberg magnet model, extending earlier results 
for the isotropic case. A large part of the paper is devoted to the 
related discrete time systems of which the pseudo-exponential poten¬ 
tial depends on an additional continuous time parameter. Our tech¬ 
niques allows us to obtain explicit solutions for the generalized discrete 
Heisenberg magnet model and evolution of the Weyl functions. 
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1 Introduction 


A skew-selfadjoint Dirac system (also called a pseudo-canonical, Zakharov- 
Shabat or AKNS system) has the form: 


d 

dx 


y(x,z) = (izj + jV(x))y(x,z) 


(x e M+, z 6 C), (1.1) 


i mi o 

0 -An 2 


V 


0 v 
v* 0 ’ 


( 1 . 2 ) 


where C stands for the complex plane, M + denotes the non-negative real 
semi-axis, I mk is the m k x identity matrix, v(x) is an m\ x m 2 matrix 
function, which is called the potential of the system, and j and V(x) are 
m x m matrices, m := m± + m 2 . Note that (jV)* = —jV, and therefore the 
system ( 1 . 1 ) is called skew-selfadjoint. 

Like a selfadjoint Dirac system fj^y(x,z) = i(zj + jV(x))y(x, z), the sys¬ 
tem ( 1 . 1 ) is also an auxiliary system for various important integrable non¬ 
linear wave equations and the case mi 7 ^ m 2 corresponds to rectangular and 
multicomponent versions of these equations. Here we solve explicitly (in 
terms of Weyl functions) direct and inverse problems for the system (1.1) for 
the case when the rectangular (possibly non-square) potential v is pseudo¬ 
exponential (see formula (2.1)). The direct problem consists in constructing 
the Weyl function and inverse problem is the problem to recover v from the 
Weyl function. 

We also derive explicit solutions of direct and inverse problems for a 
discrete analogue of the skew-selfadjoint system ( 1 . 1 ), namely, for the system 


Vk+ 1 ( 2 ) = (An + iz 1 C k ) y k {z), C k = U* k jU k , where 

UfU k = U k U* k = I m , k = 0,1,2,.... (1.3) 


This system was studied in [27] for the important subcase m x = m 2 and ex¬ 
plicit solutions of direct and inverse problems were obtained for that subcase. 
See [18] for the discrete analogue of selfadjoint Dirac system (and its general 
Weyl theory). 

A large part of the paper is devoted to the case when the system (1.3) 
depends on an additional continuous time parameter. A special choice of the 
additional parameter allows us to introduce a generalized discrete Heisenberg 
magnet model. We use our results on system (1.3) with a general-type j in 


2 






order to construct explicit solutions and the evolution of the Weyl function 
for this generalized model. The results obtained generalize those for the case 
m\ = m 2 = 1 in [27] which dealt with the well-known discrete isotropic 
Heisenberg magnet model [14,47]. 

History. Direct and inverse spectral problems for Dirac systems have a long 
and interesting history. For the selfadjoint case with a scalar continuous 
potential v it all started with the seminal paper [30] by M.G. Krein (see also 
the discussions in [3,45]). This theory is closely related to the Weyl theory 
(see, e.g., [3,8,34,39,41,44] and references therein). The Weyl theory for the 
skew-selfadjoint systems of type (1.1) was developed much later [8,16,23,35, 
40,41] and the method of operator identities [42-44] played a fundamental 
role in these studies. 

Particular subclasses of general-type potentials v are of special interest, 
since, for some subclasses, direct and inverse problems may be solved ex¬ 
plicitly. Explicit solutions of spectral problems for selfadjoint Dirac systems 
with the strictly pseudo-exponential square potentials v were given in [1,2]. 
These constructions were based on the procedure to solve direct and in¬ 
verse problems for systems with general-type potentials. A few years later 
direct and inverse problems in terms of spectral and Weyl functions were 
solved explicitly in [22] for a wider class of selfadjoint Dirac systems, that 
is, for systems with pseudo-exponential square potentials v (see also [24]). 
The same problem for system (1.1) was dealt with in [23]. Moreover, the 
above mentioned problems for systems with pseudo-exponential potentials 
were studied directly instead of using procedures for general-type selfadjoint 
(or, correspondingly, skew-selfadjoint) Dirac systems as in [1,2]. We note 
that direct methods in explicit solution of spectral problems go back to the 
works by Backlund [4] and Darboux [10,11] with essential further develop¬ 
ments in [9,12,19,20,31]. Finally, direct and inverse problems for selfadjoint 
Dirac systems with the pseudo-exponential rectangular (possibly non-square) 
potentials v were solved explicitly in terms of Weyl functions in [17]. 

Contents. The paper consists of four sections (the present introduction in¬ 
cluded) and an appendix. Section 2 deals with direct and inverse problem 
for continuous time systems. In this section we develop further the results 
from [23] on skew-selfadjoint system (1.2) with square potentials, also using 
some ideas and results from the works [16,17] on systems with general-type 
rectangular (non-square) potentials. Section 3 treats direct and inverse prob¬ 
lems for discrete time systems, generalizing earlier results from [27] to the 
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rectangular (non-square) case. Finally, Section 4 is devoted to generalized 
discrete Heisenberg magnet model, which is equivalent to the compatibility 
condition of two systems depending on two parameters t and k where t is 
non-negative real and k £ No- For the convenience of the reader the appendix 
Section A presents a number of general facts regarding admissible quadru¬ 
ples that are used throughout the paper. This theory is closely related to 
the triples approach [22-24,27] and generalized Backlund-Darboux transfor¬ 
mation (GBDT) method [15,17,37,38,41] (see also references therein). 

Notation. We conclude with some information on notations that are used 
throughout the paper. As usual N 0 stands for the set of non-negative integers, 
i.e., the natural numbers with zero included, and R + denotes the set of non¬ 
negative real values. The symbols C + and C_ stand for the upper and lower 
half-plane, respectively. Furthermore, C + denotes the closed upper half-plane 
(i.e., C_|_ U M), and Cm stands for the open half-plane {z : 5s (z) > M > 0}. 
By || • || we denote the £ 2 vector norm or the induced matrix norm, Span 
stands for the linear span, and a (a) stands for the spectrum of a. The class 
of m 2 x Wi contractive matrix functions (Schur matrix functions) on some 
domain 11 is denoted by S m2Xmi (k2). We write S > 0 when the matrix S is 
positive definite. The matrix (ct^ 1 )* is denoted by a~* and Irn a stands for 
the image of a. 

2 Continuous case: direct and inverse 
problem 

We begin with introducing the notion of a pseudo-exponential potential. The 
starting point is a matrix function v of the form 

v(x) = 2 , dle ixa *S(x)~ 1 e ixa 'd 2 , x £ K+. (2.1) 

Here i9i and $2 are matrices of sizes n x m\ and n x m 2 , respectively, a is an 
n x n matrix, and S is the n x n matrix function given by 

PX 

S(x) = S 0 + A{t)jA(t)*dt, S 0 > 0, A(x) := [e~ ixa i3 1 e ixa d 2 ] . ( 2 . 2 ) 
Jo 

Furthermore, we require the initial value So in (2.2) to satisfy the following 
matrix identity: 

a So - S 0 a* = + $ 2 $ 2 )• ( 2 -3) 
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From (2.1) and (2.2) it is clear that the potential v in (2.1) is uniquely 
determined by the quadruple {ct, S 0 , $ 1 , $ 2 }. combination with S 0 being 
positive definite, the identity (2.3) implies that S(x) is also positive definite; 
see (2.8) below. In particular, S(x) is invertible for each x > 0, and hence v 
is well-defined. 

When S 0 > 0 and (2.3) holds we call the quadruple {a, S 0 , $ 1 , $ 2 } an 
admissible quadruple. We call v in (2.1) the pseudo-exponential potential 
generated by the admissible quadruple {a, So, $ 1 , $ 2 }- See appendix Sec¬ 
tion A for a brief review of properties of an admissible quadruple and the 
relation with the theory of S-nodes. 

The definition of a pseudo-exponential potential given here is somewhat 
different from the definition in [23] which starts with [23, Eq. (0.2)]. However, 
Proposition 1.1 in [23] tells us that a pseudo-exponential potential in sense 
of [23] is also a pseudo-exponential potential as defined above. With some 
minor modifications the reverse implication is also true (see Proposition 2.8 
at the end of this section). 

In the present paper, as opposed to [23], we do not require the matrices 
$1 and $2 to be square, i.e., mi and m 2 are not required to be equal. Note 
that in [23] the matrix So is just the n x n identity matrix but, as the next 
lemma shows, it is convenient to allow So to be just positive definite. 

Lemma 2.1. Let {a, So, $ 1 , $ 2 } be an admissible quadruple, and for each 
x G M+ let E(x) be the quadruple defined by 


£(x) = {a,S(x),e- l * Q ?A,e“ a r? 2 }, 


(2.4) 


where S{x ) is the matrix defined by the first identity in (2.2). Then E(x) is 
an admissible quadruple for each x G M + . 

Proof. We first show that 


aS(x ) — S(x)a* = iA(x)A(x)*, x G M+. 


(2.5) 


To prove this identity note that the relations in (2.2) are equivalent to 


—A(x) = —iaA(x)j, A(0) = [$1 $ 2 ] (^ £ M+); 

dx L J 

■j~S(x) = A(x)jA(x)*, 5(0) = So (i£ R+). 


( 2 . 6 ) 


(2.7) 


5 


Here j is the signature matrix defined (1.2). If follows that 

— (aS(x) — S(x)a*j = aA(x)jA(x)* — A(x)jA(x)*a*, 

-^-iA(x)A(x)* = aA(x)jA(x)* — A(x)jA(x)*a*. 

On the other hand, using (2.3), the functions aS(x) — S(x)a* and iA(x)A(x)* 
have the same value at x — 0. But then (2.5) holds true. 

It remains to show that S(x) is positive definite for each x > 0. To do 
this we prove the following inequality: 

e~ ixa S{x)e ixa * > S 0 , x G M+. (2.8) 

Using (2.7) and (2.5) we have 

(e~ ixa S (x)e ixa *^J = e ~ ixa (J-S(x)\ 

- ie~ ixa (aS(x) - S(x)a*) e ixa * 

= A(x)jA(x)* + A(x)A(x)* 

= 2e~ ixa ^ 1 e ixa * >0, x G R+. 

But then 



In particular, (2.8) holds true. □ 

We note that the quadruple {a, So, $i, $ 2 } generating v coincides with 
E(0) in (2.4). 

Using straightforward modifications of the proof of [23, Theorem 1.2] (or 
particular cases of the more general [37, Theorem 3] or [38, Theorem 1.2 and 
Proposition 1.4]) we obtain the next proposition. 

Proposition 2.2. Let v be the pseudo-exponential potential generated by the 
admissible quadruple E( 0 ) = {a, So, $ 1 , $ 2 }, and let £(x) be the admissible 
quadruple defined by (2.4). Then the fundamental solution u of the system 
( 1 . 1 ), normalized by w( 0 , z ) = I m , where m = + m 2 , admits the following 

representation 

u(x,z ) = W s(x) (z)e ixzj W m (z)-\ (2.9) 

Ws(x)( 2 ) := I m + iA(x)*S(x)~ 1 (zI n - a ) -1 A(x). (2.10) 
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The transfer function of the form (2.10) (the transfer function in Lev 
Sakhnovich form) was introduced and studied in [42] (see also [41,43,44] and 
references therein). In our case, we refer to Wy,{ x ) as the the transfer function 
associated with the admissible quadruple £(x); see the first paragraph of 
Section A. The next result is a simple generalization of [23, Proposition 1.4]. 

Proposition 2.3. Pseudo-exponential potentials v are bounded on the semi¬ 
axis x > 0. 

Proof. Assume that the potential v is generated by the admissible quadruple 
£(0) = {a, So, $i, $ 2 }, and let X(x) be the admissible quadruple defined by 
(2.4). From Lemma A.l we know that the eigenvalues of a belong to C + . 
Furthermore, using the identity (A.9) with E(x) in place of X, we see that 

i(z — z)A(x)*(zI n — a*)- 1 S(x)~ l (zl n — a)^ 1 A(a;) < J m , zE C_. (2.11) 

Since the resolvent {zl n — a) -1 is well-defined in any open domain in C_, 
we have Span zen (zl n — a)^ lr dk A $k f° r k = 1,2. Therefore, inequality (2.11) 
(where A(x) is given by (2.2)) implies that 

sup ( ||S(x)- 1/2 e-^i|| + ||S(x)- 1/2 e^ 2 || ) < 00 . (2.12) 

x>0 ' ' 

From (2.1) and (2.12), it is immediate that v is bounded, that is, for some 
M > 0 we have 


n(x)|| < M x £ R + , 


(2.13) 


which completes the proof. 


□ 


Weyl function: the direct problem. The concept of a Weyl function 
of Dirac system has a long history (see the Introduction). Following the 
definition of a Weyl function for Dirac systems with square potentials (see 
also [41] for the case of non-square potentials) we say that a meromorphic 
function p satisfying (2.13) is a Weyl function of the system (1.1) whenever 
it satisfies the inequality 



u(x, z)*u(x , z) 


1 mi 

<p(z) 


dx < 00 , 


Z G C M- 


(2.14) 
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Here u is the fundamental solution of (1.1) normalized by w(0, z) = I m . We 
note that for skew-selfadjoint Dirac systems (1.1) Weyl functions have been 
introduced in [16] in an equivalent but different way. However, Proposition 
2.2 and Corollary 2.8 from [16] immediately yield the following result con¬ 
cerning direct problem. 

Proposition 2.4. Let (2.13) hold. Then there is a unique function p such 
that (2.14) is valid. This function p is analytic and contractive in Cm (he., 

T ^1712 xmi (C M )). 

For the case of a pseudo-exponential potential v we produce an explicit 
expression for the Weyl function. 

Theorem 2.5. Consider the skew-selfadjoint Dirac system (1.1), and assume 
thatv is a pseudo-exponential potential generated by the admissible quadruple 
£(0) = {a, <So, $i, $ 2 }- Then the Weyl function p of the Dirac system (1.1) 
is given by 

p(z) = — a x ) -1 , #i, a x a — ididlSQ 1 . (2.15) 


Proof. From the proof of Lemma A.2 we know that (2.15) can be rewritten 
in the form: 


ip(z) = id* 2 S 0 1 (zl n - a) 1 d 1 (l mi + id* 1 S 0 1 (zI n -a) *$ 1 ) \ (2.16) 


Taking into account the equivalence of (2.15) and (2.16) and using (2.10) 
together with the second equality in (2.2) (both at x = 0), we derive that p 
of the form (2.15) satisfies the relation: 


r 1 1 


17 1 

J m 1 

p(z)_ 

— We(o)(z) 

J mi 

0 


-i„a t 1 


(l mi + i'd* 1 S 0 1 (zI n -a) 1 di) 


Formulas (2.9) and (2.17) imply that 


u(x, z) 


1 m\ 

p(z) 


= e ixz W E(x) (z) 


L 777-1 


(2.17) 


( 4 , 1 +«;s 0 -v„-a)- 1 tfi) ‘. (2.18) 


Using Sq > 0 and taking inverses in (2.8), we see that 










This yields the inequality 

sup ||e _ ’ 7X S'(x) _1 || < oo for sufficiently large values of rj. (2.19) 

x>0 

From (2.19) and the identity (A.9) with E(x) in place of E, it then follows 
that the function x —> e^ z ~ z ^ x Ws( x )(z)*Ws^ x j(z) belongs to L^ xm (0, oo) for 
sufficiently large values of A(z): 

6 L 2 mxm ( 0, oo). (2.20) 

That is, the entries of & i( '' z ~ z ' )x Wy j ( x )[z)*Wt,(x)(z) are squarely summable with 
respect to x for sufficiently large values of £y(z). Finally, in view of (2.18) and 
(2.20), the inequality (2.14) holds for <p given by (2.15) and for sufficiently 
large values of A(z). 

Taking into account Proposition 2.4 and the analiticity of ip given by 
(2.15), we see that (2.14) holds for all z G (Cm \ cr(a: x )), that is, ip is the 
Weyl function. □ 

Weyl function: the inverse problem. Theorem 2.5 presents the solution 
for the direct problem. We now turn to the inverse problem. The uniqueness 
theorem below is immediate from [41, Theorem 3.21 and Corollary 3.251 (see 
also [16]). 

Theorem 2.6. Let ip be a Weyl function of a skew-self adjoint Dirac system 
with a potential v which is bounded on [0, oo). Then this v can be uniquely 
recovered from ip. 

For the case of pseudo-exponential potentials we have an explicit proce¬ 
dure to recover the potential from the Weyl function. This procedure uses 
such well-known notions from control theory as realization, minimal realiza¬ 
tion and McMillan degree (see, e.g., [5,28] or [41, App. B]). We note that, 
according to Theorem 2.5, the Weyl function ip of a skew-selfadjoint Dirac 
system with a pseudo-exponential potential is a strictly proper rational ma¬ 
trix function. 

Theorem 2.7. Let p be a strictly proper rational m 2 x m\ matrix func¬ 
tion. Then ip is the Weyl function of a skew-selfadjoint Dirac system with 
a pseudo-exponential potential v generated by an admissible quadruple. The 
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corresponding quadruple {a, So, t?i, $ 2 } can be obtained explicitly by using the 
following procedure. First, we construct a minimal realization of p : 

p{z)=i0l{zl n - 1 )- l e l . (2.21) 

Next, we choose X to be the unique positive definite solution X of the Riccati 
equation 7 X — X 7 * — iX9 2 9%X + i9\9* = 0. Finally, we put 

S 0 = I n , = X~ l,2 9i, d 2 = X 1 ' 2 9 2 , a = X-V^X^ + Uhtil. (2.22) 

Proof. Given a strictly proper rational matrix function p, Theorem A.3 
provides the described above three step procedure to construct a quadru¬ 
ple {a, So, $ 1 , $ 2 } such that representation (2.15) of ip holds. It also fol¬ 
lows from Theorem A.3 that the procedure is well-defined and the quadruple 
{a, So, $ 1 , $ 2 } is admissible. Then we know from Theorem 2.5 that p we 
started with is precisely the Weyl function of the skew-selfadjoint Dirac sys¬ 
tem with the potential v generated by {a, So, d 1 , $ 2 }- □ 

As we note in the second paragraph before Lemma 2.1 the definition of 
a pseudo-exponential potential given in the beginning of this section differs 
from the definition employed in [23] which starts from formula (0.2) in [23]. 
The next proposition presents the analogue of formula (0.2) in [23], which 
coincides with (0.2) for the case mi = m 2 and So = I n . The conclusion is 
that the two definitions lead to the same class of potentials. 


Proposition 2.8. Let v be the pseudo-exponential potential generated by the 
admissible quadruple {a, So, $ 1 , $ 2 }, and let A be the 2 n x 2 n matrix defined 
by 


a* 0 
a 


(2.23) 


Then the potential v is also given by 


v{x) = 2 $*( [So 


-Un] e~ 2ixA 




x G M+. 


(2.24) 


The proof of the proposition is given in Appendix (and is close to the 
considerations in [23]). 
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3 Discrete case: direct and inverse problem 

Recall that the discrete skew-self adjoint Dirac system (SkDDS) has the form: 

y k + 10 ) = (I m + iz~ x C k ) y k (z ), C k = U k jU k (k G N 0 ). (3.1) 

Here No stands for the set of non-negative integers, the matrices U k are 
unitary of size m x m, where m does not depent on k, and j is the m x m 
signature matrix defined by ( 1 . 2 ) with m = + m 2 and with m\ and m 2 

not depending on k. The sequence is called the potential of the 

system. Note that the second part of (3.1) implies that C k = C k = C k l for 
each k G N 0 . 


Definition 3.1. The Weyl function of SkDDS is an m\ x m 2 meromorphic 
matrix function <p(z) on Cm (for some M > 0), which satisfies the inequality 


OO 

E W( Z Y W k (z)*W k {z ) 

k=0 


<P(z) 

Im?. 


< OO, 


(3.2) 


where w k (z) is the fundamental solution of SkDDS normalized by w 0 (z) = I m . 

We shall consider the case when the potentials {Cfc}fc e N 0 are pseudo¬ 
exponential (see Definition 3.3 below), and we shall show that for such po¬ 
tentials the Weyl function can be constructed explicitly. 

Pseudo-exponential potentials have been introduced in [27] for the square 
case when m\ = m 2 . We shall show that the same scheme developed in [27] 
for constructing the corresponding Weyl function also works in the non¬ 
square case when m\ and m 2 are not equal. 

Similar to the continuous case, our starting point to define pseudo-expo¬ 
nential potentials is an admissible quadruple S 0 = {cc, Ro; ^15 ^ 2 }, and A 0 = 
[$1 t^]. I 11 addition, we assume that a is non-singular. We set 

Afc+i = A k + ioL~ l k k j (k e N 0 ); (3.3) 

S k+ 1 = S k + a~ l S k a~* + a~ l K k jK* k a~* (.k e N 0 ), (3.4) 


where a * stands for (a x )*. From (3.3) and A 0 = [ih 


■d 2 ] it follows that 


Afc = [(I n + ia~ l ) k Di (I n - ia- 1 )^] (k G N 0 ). (3.5) 
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By definition E fc is the quadruple given by 

= {a, S k , ( I n + ia~ l ) k d i, (I n - ia~ l ) k $ 2 }- (3.6) 

The next lemma shows that this quadruple is again admissible if, in addition, 
the pair {a, $ 1 } is controllable. 


Lemma 3.2. Let So = {a, So, $ 1 , $ 2 } be an admissible quadruple. Assume 
that the pair {a, id 1 } is controllable. Then a is non-singular, and the quadru¬ 
ple E fc defined by (3.6) is admissible for each k G N 0 . Moreover the pair 
{a, (I n + ia~ l ) k di} is controllable. 

Proof. The fact that the pair { 0 ,$!} is controllable implies that the pair 
{a,Ao} is also controllable. But then item (ii) in Lemma A.l tells us that 
cr(a) C C + . In particular, a is non-singular, and thus the quadruple E& is 
well-defined. 

Next we use (3.3) and (3.4) to prove that 

aS k — S k a* = iA k A* k (k G N 0 ). (3.7) 

This will be done by induction. Since A 0 = [ 7 ? 1 i d 2 \ and the quadruple 
E 0 = {a, S 0 , ^ 1 ,^ 2 } is admissible, the identity (3.7) holds for k — 0 . Suppose 
now that (3.7) holds for k = r. Then, using the expression for S r+ ± from 
(3.4), we obtain 

aS r+ 1 — S r+ ia* = iA r A* + ia~ l A r A*a~*+ 

+ A r jA*a~* - a _ 1 A r jA*. (3.8) 

On the other hand, using (3.3) for k = r, we obtain 

iA r+ iA* +1 = i(A r + ia~ 1 A r j)(A* — ijA*a~*) 

= iA r A* - ct _ 1 A r jA* + A r jA* r a~* + ia _ 1 A r A*a~*. (3.9) 

Together (3.8) and (3.9) yield (3.7) for k = r +1, and thus (3.7) holds for all 
k e No- 

Next we show that S k is positive definite for each k e N 0 . To do this we 
first use that cr(a) C C+. The latter inclusion and the identities (3.7) imply 
(cf., [26] or Theorem 1.4.1 in [21]; also (A.7)) that 

1 r°° 

Sk = ^ (A I n - a)~ 1 A fc A*(A/ n - a^dX, k G N 0 . (3.10) 


12 


Now, the inequality S k > 0 is proved by contradiction. Assuming that S k 'f- 0, 
we derive (from (3.10)) the existence of a vector g ^ 0 such that g(XI n ~ 
a) L A k = 0. In particular, using the identitiy (3.5), we see that g{XI n — 
= 0, where g = g(I n + ia" 1 )^. The fact that a(a) C C + implies that 
—i is not an eigenvalue of a, and hence the matrix I n AioT 1 is non-singular. 
But then g ^ 0 because g ^ 0. On the other hand, from the Taylor expansion 
of g(XI n — a)~ 1, $i we obtain that 

goPd i = 0 (g 0 ), for all p > 0 , 

which contradicts the assumption that the pair {a, $ 1 } is controllable. Thus 
Sk is positive definite, and hence E*, is admissible. 

It remains to prove the pair {a, (J n + ia _ 1 ) fc, d 1 } is controllable. To do 
this note that 

n— 1 7i—l 

Span Im a v (I n + ia~ l ) k d\ = (I n + ia~ l ) k ^Span ImaT'di)- (3.11) 

v =0 v =0 

Since the pair {a, ??!} is controllable, the space Span (J"Zq Im oAdi is equal 
to C n . It follows that the space in the right-hand side of (3.11) is also equal 
to C n . But then the same holds true for the space in the left-hand side of 
(3.11), which implies that {a, (I n + ia~ 1 ) k, di} is controllable. □ 

Note that the conditions in Lemma 3.2 imply that the matrix Sk is non¬ 
singular for each k. This allows us to define the following sequence of matri¬ 
ces: 

C k = j + A*SV'A fc - A* k+1 S^A k+1 , k = 0,1,2,.... (3.12) 

Definition 3.3. We call the quadruple E 0 = {a, S 0 , i?i, $ 2 } strongly admis¬ 
sible if it is admissible and the pair {a,$i} is controllable. In this case we 
refer to the sequence of matrices {CfclfcgNo * n (3-12) as the pseudo-exponential 
potential generated by E 0 . 

Later we shall see (Proposition 3.6 below) that the second part of (3.1) 
is fulfilled for any pseudo-exponential potential {C k } kG fq 0 . 


Remark 3.4. We note that for any continuous pseudo-exponential potential 
v (see its definition at the beginning of Section 2) there is a strongly admissible 


13 


quadruple which generates it. Indeed, recall that in view of Proposition 2.3 
and Theorem 2.6 there is a unique solution of the inverse problem considered 
in Theorem 2. 7. Thus, starting from the Weyl function of a system with the 
pseudo-exponential v we recover some admissible quadruple generating this v 
via formula (2.22) from Theorem 2.7. It is immediate from (2.22) that the 
corresponding pair {a, $ 1 } is controllable, and so the recovered quadruple is 
strongly admissible. 

Next, we present a discrete analog of formula (2.9) for the fundamental 
solution. See the next theorem which is a minor generalization of Theorem 
0.1 in [27], 


Theorem 3.5. Let £ 0 = {a, S 0 , $ 1 , $ 2 } be a strongly admissible quadruple, 
and let {Cfej/cgN,, be the pseudo-exponential potential generated by S 0 . Then 
the fundamental solution w k (z) (k G No) of the normalized initial value prob¬ 
lem 


w k+ i(z) 


ylm + Z Ck j W k( Z )i W o(z) = 4 


(3.13) 


can be represented in the form 


Wk(z ) 


W Sk (-z) 


I'm T j 

Z 


Wx 0 (-z)-\ 


(3.14) 


where is the admissible quadruple defined by (3.6) and Ws k is the transfer 
function associated with T k , that is, 

W^ k {\) = I m + tAlS^\XI n -a)- 1 A k . (3.15) 


See the first paragraph of Section A for the definition of the transfer function 
associated with an admissible quadruple. 

Proof. The proof is based on the following identity: 


WWA) (l m - id = (l m - led W^(A). (3.16) 

This identity is the analogue of formula (2.2) in [27], and its proof follows 
the same line of reasoning as the proof of formula (2.2) in [27]. We omit the 
details. 
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Next, using (3.16), we complete the proof by induction. First notice that 
the equality (3.14) holds for k = 0. Next, assume (3.14) is proved for k = r. 
Then, using (3.13), (3.14) with k = r, and (3.16) with k = r and A = —z, we 
see that 


W r -\-\ 


- C, 


w r (z ) 


(im + % -c)j W^{-Z) (l m + ~Aj W^-Z ^ 1 

Wv r+ 1 (-Z) flm + -A (im + A) W *> (--)' 1 
\ r+1 

i) 



Thus (3.14) holds for k — r + 1, and so it holds for all k G No- □ 

Our next proposition shows that the second equality in (3.1) is valid. 


Proposition 3.6. Let E 0 = {a, So, $i, $ 2 } ^ e a strongly admissible quadru¬ 
ple, and let {C k } k ^ 0 be the pseudo-exponential potential generated by S Q . 
Then we have 


C k = U* k jU k , where U* k U k = I m (k e N 0 ). (3.17) 

Proof. According to Lemma 3.2, E fc is an admissible quadruple, and so we 
may substitute T, k (instead of £) into formula (A. 8 ) from Lemma A.l. Then 
(A. 8 ) implies that 

Wx k (\)Wz k (X)* = Im NN 0 ), (3.18) 

where A stands for complex conjugate for A. Note that (3.18) makes sense 
only when A and A are no poles of W Efc . 

It is immediate from (3.12) that C k = C k . From (3.16) and (3.18) we see 
that 

(/ m — iX 1 C k ) (Im + *A 1 C k j = A 2 (A^ + 1 )I m , A el, 

that is, C'l = I m . Thus, the equality C k = Cy l is valid, and we obtain the 
representation 

C k = U k j k U* k , (3.19) 
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where U k are unitary matrices and j k are diagonal matrices, the entries of 
which take the values ± 1 . 

Next we show that we may choose jk = j■ From Lemma A.l we know 
that that a(a) C C+. In particular, —i does not belong to a (a). We first deal 
with the case when the same holds true for i. In that case (3.18) makes sense 
for both A = i and A = — i. Let us partition the transfer matrix function 
Wz k (A) into the two blocks, the first consisting of the first mi columns and 
the second of the remaining m 2 columns: 

w Sl (A) = [(WV.M), (WV,(A)) 2 ]. 

Now, take A — —i in (3.16) and multiply the resulting identity from the right 
by W-£ k {i)*. Then, using (3.18) with A = — i, we obtain 

2(W'e, + ,H)) 1 (W'e.(*)); = Im + C,. (3.20) 

Repeating the argument with A = i in place of A = — i we see that 

= I m - c*. (3.21) 


Formulas (3.20) and (3.21) imply that 

rank(J m + C k ) < mi, rank (I m - C k ) < m 2 . (3.22) 

The existence of the representation (3.17) now easily follows from (3.19) and 
(3.22). Thus the proposition is proved when i jL cr (ck). 

It remains to consider the case when i belongs to a (a). In this case 
we approximate the original quadruple E 0 = {a, S 0 , $ 1 , $ 2 } by a set of new 
quadruples £(e) = {ck + £/ n , So, $ 1 , $ 2 }, with £ = e > 0 and £ sufficiently 
small. These new quadruples are again strongly admissible. Thus they satisfy 
the conditions of Theorem 3.5 and the additional condition 

i jL ct(ck + £l n ). 

Therefore, applying the result of the previous paragraph, (3.17) is valid also 
in this case. Taking limits for £ {, 0 one obtains (3.17) for the original 
quadruple £ 0 . □ 
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Remark 3.7. It is not difficult to see that Theorem 3.5, the identities (3.7) 
and the representations in (3.17) hold under weaker conditions. Indeed, as¬ 
sume that a is non-singular, that So = Sq, that (3.7) holds for k = 0, and 
that Afc and Sk are given by (3.3) and (3.4), respectively. Then the identi¬ 
ties (3.7) are valid for fc 6 No. Assuming additionally that det S r 0 for 
0 < r < N, one can show that the fundamental solution Wk(z) (0 < k < N ) 
of the discrete system (3.13) can be represented in the form (3.14) and that 
(3.17) is also valid for 0 < k < N — 1. 


Weyl function: the direct problem. The following result is a discrete 
analogue of Theorem 2.5. The special case m\ = m 2 of this result is a 
somewhat stronger version of Theorem 0.4 in [27]. 

Theorem 3.8. Let E 0 = {a, Sq, Di, $ 2 } be a strongly admissible quadruple, 
and let {Ufc}fc G N 0 be the pseudo-exponential sequence generated by E 0 . Then 
the discrete skew-self adjoint Dirac system (3.1) has a unique Weyl function 
(p which is given by the formula 

<p(z) = -i^Sg^zIn + / 8)~ 1 i9 2 , f3 := a - (3.23) 

Moreover, satisfies (3.2) in the half-plane Sx: > \ (a finite number of points 
excluded). 

Proof. Part 1. In order to show that is a Weyl function one first 
compares (3.23) with (A.4) to see that is the (second) function associated 
with E 0 . Thus, the second part in (A. 11 ) implies that 

cp( z ) = b(-z)d(-zy 1 , (3.24) 


where b and d are the blocks of W% 0 (see (A.10)). Using (3.24) we have 


W£o (~Z) 


-1 


if{z) 


L m 2 


= W^(-z 


l-l 


X~z) 

d(—z) 


d(-z)- 1 = 


0 


L m,2 


d(-z)- 1 . 

(3.25) 


Next, for k G No, let Ws k be the transfer function of the quadruple E*, defined 
by (3.6). Recall that the normalized fundamental solution Wk of (3.1) is given 
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by (3.14). Using (3.25) this yields 


w k [z) 


'<P(z) 


L rri 2 


= Wx t (-z)[I m +-j) W^(-z) 


-1 


'(p(z) 


L m 2 


/ • \ k 

( 1 \ 

' 0 ' 

‘ c-*, 

i ' 

+ 

g 

1 i 

Jrri2_ 

V z J 


Z — l 


■ \ k 


Wz,J-z) 


d(—z) 1 
d{—z)~ l . 


(3.26) 


Since S7 > 0 and cr(a) C C + (according to item (i) in Lemma 3.2), the 
identity (A.9) implies that 


Wv„(-zywv k (-z) < Im (z G C+). (3.27) 

From (3.26) and (3.27) we easily derive (3.2) for all z G Ci/ 2 , a finite number 
of eigenvalues excluded. 

Part 2. It remains to prove the uniqueness of the Weyl function. We shall 
again need the identity (A.9) (for the quadruple £ fc ): 

W Sfc (A)*W Sfe (A) = I m - *(A - A)A*(A I n - aT^k^XIn ~ «)“%. (3.28) 

Since the left upper block of (3.28) is nonnegative, taking into account 
(3.5), we obtain for z = (—A) G C + the inequality 

dl(zl n + a*)~\l n - i{ct*)-') k S k \l n + ia-^izln + a)” 1 ^ < —^=/ TOl . 

z z 

(3.29) 

Using the same arguments as at the end of the proof of Lemma 3.2, we see 
that the controllability of the pair {a, $ 1 } implies the equality 

Span zeOe (z 0 ){zIn + a)~ 1 '$i = C n (3.30) 

for any ^-neighborhood O £ (z 0 ) of any z 0 G C. From (3.29) and (3.30) we see 
that for all k G N 0 the positive definite matrix (J n — i{a*)~ l ) k S^ 1 {I n +ia~ 1 ) k 
is uniformly bounded. Hence, the matrix function 

i(z - z)$l(zl n + a*)~ 1 (I n - i(a*)~ 1 ) k S / 7 1 (I n + ia~ 1 ) k (zl n + a) _1, di 
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is sufficiently small for all z G Cm and sufficiently large values of M. Recall 
that the matrix function above coincides with the left upper block of the 
second term on the right hand side of (3.28) (for A = —z). Therefore, 
relations (3.14) and (3.28) yield the inequality 



0] Wx 0 (-z)*w k (z)*w k (z)Wx 0 (-z ) 



> 



(3.31) 


for z G Cm and sufficiently large values of M. Thus, for sufficiently large M 
and all z G Cm, we have 


dim!/ < m 2 , L := {g G C m | g*w k (z)*w k (z)g < oo}. (3.32) 


k =0 


Assume now that for some z G Cm there is a value ip(z) such that 


|V(*)* 7 ma] 


w k {z w k lz 


k =0 


if(z) 


L m 2 


< oo, ip(z) 7 ^ tp(z). (3.33) 


Then (3.2) and (3.33) contradict the inequality (3.32). Hence, there is no 
such tp and the Weyl function is unique. □ 


Weyl function: the inverse problem. From Theorem 3.8 we know that 
the Weyl function of a discrete skew-selfadjoint Dirac system (3.1) with a 
pseudo-exponential potential generated by a strongly admissible quadruple 
is a strictly proper rational matrix function. The next theorem shows that 
the converse is also true. 

Theorem 3.9. Let <p be a strictly proper rational mi x m 2 matrix func¬ 
tion. Then p is the Weyl function of a discrete skew-selfadjoint Dirac system 
(3.1) with a pseudo-exponential potential generated by a strongly admissible 
quadruple S 0 . A corresponding S 0 = {a, So, $ 1 ,^ 2 } can be obtained explicitly 
by using the following procedure. First, we construct a minimal realization 

of V ■ 

(p(z) = —id*(zl n + 7 )~ 1 d 2 ■ (3.34) 

Next, we choose X to be the unique positive definite solution of the Riccati 
equation jX — Xj* — iXOid^X + id 2 0\ = 0. Finally, we put 

So = In, 0 1 = X^di, d 2 = X- l ' 2 d 2 , a = X-V^X 1 ' 2 + 10201. (3.35) 
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Proof. Given a strictly proper rational matrix function ip, Corollary A.4 
shows that the procedure to recover E 0 = {a, S 0 , $ 1 , r? 2 } is well-defined and 
So is strongly admissible. Moreover, according to Corollary A.4, p admits 
representation (3.23) where the quadruple {a, So, $i, $ 2 } i s gi ven by (3.35). 
Then we know from Theorem 3.8 that p is the Weyl function of system (3.1) 
with the potential { C k } generated by {a, S 0 , i?i, $ 2 }- □ 

We conclude this section with some auxiliary results on the mxm matrices 
H k and H k \ 


nt :=2W^(i)P I W Sl! (-,r, P, = (I m +j )/ 2 ; (3.36) 

K ~ 2W Eh (-i)P 2 W Ei ( i)\ P 2 = ~ j)/ 2, (3.37) 

which will be essential in the next section. Here Ws k is the transfer function 
of the form (3.15). 

Lemma 3.10. Let E 0 = {a, S 0 , $ 1 , $ 2 } be a strongly admissible quadruple, 
let E fc be defined by (3.6), assume that i ^ &(a), and let {Cfc}fc e N 0 be the 
pseudo-exponential potential generated by E 0 . Then 

(Im + Cfc) H k = H k+1 (I m + C'fe) — 0, (3.38) 

(I m ~ C k ) Hi = Ht +1 (I m - C k ) = 0. (3.39) 

Proof. Since E 0 = {ct, S 0 , $ 1 , $2} is admissible, we know (see Lemma A.l) 
that —i ^ a (a). By assumption i ^ a ( ck ). This allows us to apply (3.18), 
first with A = i and next with A = —1. It follows that both MAy (i) and 


W-£ k (—i) are well-defined and invertible. Moreover, we have 

W^(i)-' = W^(-ir and wyt-i )- 1 = W Et (i)\ (3.40) 

Using these identities the formulas for H, and can be rewritten as: 

Hi = 2 W El ( ! )c 1 ny(*)-\ n- k = 2W Ek (-i)P 2 W Ek (-i)- 1 . (3.41) 

Next, applying (3.16), first with A — i and next with A — —i, we see that 

(Im - cyny(i) = 2Wy +1 (i)P 2 , (3.42) 

(Im + C t )W Eh (-i) = 2W Eh P-i)P,. (3.43) 
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It follows that 


(4, - C k )H+ = 2 (I m - C»)W' & (i)P,W 5 ; ,(-i)* 

= 4 lV El+1 (i)P 2 f>,lV E ,(-i)* = 0, (3.44) 

(4, + C k )Hy = 2(I m + C k )W Ek (-i)P 2 W Et (i)- 

= 4lV El+1 (-i)P 1 f 2 lV El (i)- = 0. (3.45) 

From (3.41), we obtain 

ntW Ek (i) = 2W Ek (i)P 1 W Et (i)~ 1 W Ek (i) = 2 W Ek (i)P h (3.46) 

U- k W Ek (-i) = 2W Ek (-i)P,W E „(-i)- 1 W Ek (-i) = 2 W Ek (-i)P 2 . (3.47) 
In particular, we have 

'HtW^ k (i)P 2 = 0 and H^Wx k (-i)Pi = 0 (k > 0). (3.48) 

Using (3.42) and (3.43) again and taking into account (3.48), we see that 

Ht +1 (I m - C t )W Ek (i ) = H* +1 (I m - C t )W Ek (t)(P t + P 2 ) 

= 2 nt +k w Ekkt (i)p 2 (p 1 + p 2 ) 

= 2Ht +1 W Ekkt (i)P 2 = 0, (3.49) 

and 

H~ k+1 (I m + C k )W Ek (-i) = H~ k+1 (I m + C k )W E _„(-i)(P, + P 2 ) 

= 2 ny +1 w Ek+t (-i)p 1 (p 2 + p 1 ) 

= m- k+1 W EM (-i)P k = 0. (3.50) 

Formulas (3.44), (3.45), (3.49) and (3.50) yield (3.38) and (3.39). □ 

For later purposes (see the next section) we mention that under the as¬ 
sumptions of Lemma 3.10 the matrices and can be rewritten as 

ut = i m + w Ek (i)jw Ek (-iy, (3.5i) 

HI =I m -W Ek (-i)jW Ek (i)‘. (3.52) 

Indeed, these identities follow from (3.36) and (3.37) by using the two inden- 
ties in (3.40). Together (3.51) and (3.52) imply that 

Hi + (H- k y = 24 ,. 
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Finally, using (3.51) we see that 


Ht = I m + (I m + iA* k S k \iI n - «)“%) 3 

x (l m + zA* k S^ 1 (-zI n -a)- 1 A k y 
— Im + (j + iA-lS k 1 (il n — ot) 1 A fc j) 

x [l m - iA* k {iI n - a*)~ 1 S k 1 A k ) 

= Im+j + iAIS^iiln - a)~ 1 A k j 

- ijA* k (H n - a*)~' 1 S k 1 A k 

+ A* k S k \iI n - a)- l A k ]A* k (tI n - a*)- 1 S^ 1 K k . (3.53) 

Analogously, using (3.52), we have 

U k =Im- (Im + lAlS k \-lI n ~ a)^ A k ) J 

x (I m + iAlS^iiln - a)~ 1 A k )* 

= Im - {j + iA* k S k l (-iI n - a) _1 A fc j) 

x (I m - iA* k (-iI n - a*)~ 1 S k 1 A k ) 

= I m - j + iAlS^iiln + a)~ l A k j 
-ijAKiy + aT'S^Ak 

- A* k S k l {iI n + a)~ 1 A k jA* k (iI n + a*) -1 ^ (3.54) 

4 Generalized discrete Heisenberg magnet 
model and its auxiliary linear systems 

This section is devoted to the generalized discrete Heisenberg magnet model. 
We apply the quadruples considered in the previous section in order to con¬ 
struct explicit solutions of this model. Given potential {C k (t)} from such 
solution and using the results from the previous section, we also express ex¬ 
plicitly evolution of the Weyl function of the auxiliary linear system, which 
coincides with system (3.1). 

4.1 Generalized discrete Heisenberg magnet model 

For the case that mi = m 2 = 1, system (3.1) is an auxiliary linear system 
for the integrable isotropic Heisenberg magnet model [47] (see also [14, Part 
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II, Section 1.2]). More precisely, it was shown in [14,47] that the isotropic 
Heisenberg magnet model is equivalent to the compatibility condition 



Fk+iGk ~ G k F k 


of the auxiliary systems 


Uk+i — GkVk, Gk(t,z) 
-jlUk = F k y kl F k (t,z) 


Im + ~C k (t); 

Hj(t) Hi(t) 
z+i z — i’ 


where m = 2, 


Hf(t) = f k (t)(I 2 ± C k (t))(I 2 ± C k -!(t)), 


(4.1) 


(4.2) 

(4.3) 


(4.4) 


f k is a scalar function depending on C k _\ and C kl k > 0, C k has the form 


C k = 


r 3 

c k 


•2 


Cl + 1C 


4 ~ icl 

-r 3 

Gh. 


(cl) 2 + (cl) 2 + (cl) 2 = 1, (4.5) 


and the numbers c k , cl and cl are real valued. It is easy to see that rep¬ 
resentation (4.5) is equivalent to the conditions C k = C k , detCT, = 1 and 
Tr C k = 0, where Tr means trace. These conditions are in turn equivalent to 
the conditions C k = Cl and the eigenvalues of C k equal 1 and —1. Hence, 
(4.5) is equivalent to the second part of (3.1), where mi — m 2 — 1 (see [27]). 
Moreover, the second part of (3.1) yields relations 


I 2 + C k = 2U k P 1 U* k , I 2 — C k 


2U k P 2 Ul, 


(4.6) 


where Pi = (I 2 + j)/2 , P 2 = (J 2 — j)/2 and P\P 2 = P 2 P\ = 0. Therefore, 
(4.4) and (4.6) imply that H k has rank 1, that 

(h TC k (t))H*(t ) = Hi(t){h =F C k -,(t)) = 0, (4.7) 

and that (4.7) defines matrices H k and H k up to scalar factors. 

Generalized Heisenberg magnet models are also actively studied (see, 
e.g., [13,25,29,46,48,49]), including the case of the continuous generalized 
Heisenberg magnet model, where the mxm matrix function C (an analog of 
C k ) satisfies the condition C = C* = C~ x (see [13, Remark 5] and references 
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therein). Using an analog of condition (4.7) instead of (4.4), we consider the 
following generalized discrete Heisenberg magnet model (DGHM): 


ijC k (t) = (H k+ ,(t) - H* +1 (t)) C„(t) 

-C k {t)(Hi(t)-H+{t)), (4.8) 

where 

(4 + C»(t)) H k (t) = H k+1 (t) (4 + C k (t)) = 0. (4.9) 

(4 - C k (t)) H+(t) = ff+ +1 («) (4 - C k (t)) = 0, (4.10) 

and 

c k (t) = U k (t)jU k (t)*, U k (t)U k (ty = I m , k G No. (4.11) 

The next proposition shows that (4.8) coincides with (4.1). 

Proposition 4.1. Under conditions (4.9)-(4.11), the discrete generalized 
Heisenberg magnet model (4.8) is equivalent to the compatibility condition 
(4.1) of the auxiliary systems (4.2) and (4.3). 

Proof. In view of the definitions of G k and F k , we rewrite (4.1) in the form 

Im + C'k 

Z 

(4.12) 

Using the equality z (t±^ — h ~ after simple transformations we see that 
(4.12) is equivalent to the relation 

4 C * ~ H: +1 )c k - C k (H- k - H+) 

+ ^ ^ {(Im — C k )Ht — H+ +1 (I m — C k )) 

+ - A((/ m + C k )H~ - H~ +1 (I m + C k )). (4.13) 

Z X 

Taking into account the equalities (4.9) and (4.11) , we derive that the terms 
in the second and third lines of (4.13) both turn to zero, and therefore (4.13) 
is equivalent to (4.8). □ 



~Ck= (lm + -C k 
z at \ z 


Hi 


z + i 


- + 


Hi 


II 


fc+i 


H. 


z + i 


- + 


fc+i 
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In order to construct explicit solutions of DGHM, we introduce an addi¬ 
tional parameter t into the quadruples generating potentials (Cfc(f)}. Namely, 
we consider quadruples E f = E 4)0 = {a, So(t), Here a is con¬ 

stant and the dependence of the quadruples on t is described by the equations 

(W) (t)=-2((a-il n ) 1 A 0 (t)Pi + (a + iI n ) 1 A 0 (t)P2), (4-14) 

Pi : = (4n + j)/ 2, P2 ■= (Im - j)/ 2; (4.15) 

^ = ~ + ( a + 

+ So(t)(oc + il n ) 1 + S'o(t)(cr — il n ) 

+ 2(a 2 + I n ) 1 (aAo(t)jAo(t)* + Ao(t)jAo(t)*a7) 

xdo^ + y- 1 ). (4.16) 

Formula (4.14) yields the equalities 

A 0 (t) = [i?i (t) $2 (t)] , (4.17) 

t?i(t) = e- 2 t(Q -^ )_ Vi(0), tf 2 (t) = e- 2 ' (Q+i 7 ™ ) ' 1 ^ 2 (0). (4.18) 

We assume for simplicity that E 0 ,o H strongly admissible and that i qL er(a). 
Then, the quadruples E 4j0 are well-defined and uniquely determined by the 
initial quadruple E 0i o- Moreover, all E t>0 are strongly admissible (see Subsec¬ 
tion 4.2). It was shown in Section 3 that the strongly admissible quadruples 
S f 0 determine strongly admissible quadruples E <>fc using relations (3.3) and 
(3.4), and that in this way E t)0 generate potentials (C'fc(t)}. Since E t)0 are 
determined by E 0 ,o we also say that E t ,k are determined by E 0 ,o and that 
potentials {Ck{t)} are generated by Eo,o- The following theorem is our main 
result in this section. 

Theorem 4.2. Let Eo,o = {a, So, $ 1 , $ 2 } be a strongly admissible quadruple, 
and assume that i a (a). Then for all t e K and k E N 0 the quadruples 
E t) fc determined by E 0i o are strongly admissible. Furthermore, the matrix 
functions 

Ht(t) :=«+(*) = Wx v (i)(I m + j)Wz v (-iy, (4.19) 

H;(t) = Wz lf (-i)(I m - (4.20) 

are well-defined and satisfy DGHM (4.8) and conditions (4.9) and (4.10). 
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Recall that matrix functions TLf: were introduced and studied in the last 
part of Section 3. The proof of Theorem 4.2 and some related results on 
discrete Dirac systems depending on an additional continuous time parameter 
are given in the next subsection. 

Finally, we note that in view of (4.18) we have the following corollary of 
Theorems 3.8 and 4.2. 

Corollary 4.3. Under the conditions of Theorem f.2, the evolution of the 
Weyl function (pft, z ) of the system (4.2) is given by the formula 

<p(t,z) = S 0 (t)-\zl n + Z^))-^- 2 ^^ -1 ^, 

P{t) := a - j e - 2 *(«+i/„)- 1 7?2l? * e - 2 t(a*-iJ„ ) - 1 iS -o( t )- 1 j dk 0 fc ( O ) = 1>2 ). 

4.2 Discrete Dirac systems depending on an additional 
continuous time parameter 

Let £ = {a, So, a s t ron gly admissible quadruple, and assume that 

i (f a(a). We shall associate with £ a family of strongly admissible quadru¬ 
ples depending on two parameters t and k, where t E M+ and k E No- This 
will be done in two steps. 

Step 1: CONSTRUCTION of £ t . In this first step we just require £ to 
be admissible, not necessarily strongly admissible. From item (i) in Lemma 
A.l we know that —i (jL cr(a). Together with the assumption that i (jL cr(a) 
this implies that both a — il n and a + il n are invertible. This allows us to 
define A(t) = A 0 (t) and S(t) = Soft) using (4.18) (and (4.17)) and (4.16), 
respectively. 

Lemma 4.4. Let £ = {a, So, 0i, 02} be an admissible quadruple, and assume 
that i a (a). Let fijft), j = 1,2, and S(t) be defined by (4.18) and (4.16), 
respectively. Then on some interval interval —£\ <t < e 2 the quadruple 

£ t = {a,S(t),0!(t),0 2 (t)} (4.21) 

is admissible too. Moreover, ifY, is strongly admissible, then the same holds 
true for the quadruple £* for each t E M. 

Proof. Since S 0 is positive definite, the same holds true for Sft) provided t 
is sufficiently small. Put 

T (t) := a Sft) - Sft) a* - iA(t)A{t)*. (4.22) 


26 


Here A (t) is defined by (4.17). Since {a, S 0 , $ 1 , $ 2 } is an admissible quadru¬ 
ple, it follows that Y(0) = 0. After some easy transformations, using (4.17) 
and (4.16) and differentiating the right-hand side of (4.22), we see that Y 
satisfies the linear differential equation 

jT(t) = -(« + i/ n )- 1 (aY(f) + Y(t)a*)(a* - il n )~ l + 

-{a- iln )- 1 (aY (t) + Y(t)cU) {a* + il n )~\ (4.23) 

Since the initial value Y(0) is zero, it follows that Y(t) = 0, that is, 

aS{t) - S(t)a* = iA(t)A(t)*. (4.24) 

Using (4.18) and (4.17), we conclude that {a, S(t), $i(t), ^(i)} is an admis¬ 
sible quadruple. 

Next, assume £ is strongly admissible. The identity (4.24) was already 
proved above. By Definition 3.3, the fact that £ is strongly admissible means 
that the pair a and d\ is controllable. But then the first equality in (4.18) 
tells us that 

n— 1 Ti—l 

Span |J Ima^i(f) = e-^D-^U 1 (^pan (J i m (^(O)) = C, (4.25) 

v=0 u=0 

that is, the pair a and 'di(t) is also controllable. Furthermore, since £ is 
strongly admissible item (ii) of Lemma A.l shows that cr(a) C C+. From 
(4.24) and o’(a) C C + follows representation (A.7) of So(t) = S(t). Since the 
pair a and $i(t) is controllable representation (A.7) implies that S(t) > 0. 
Hence, again using Definition 3.3, the quadruple £/ is strongly admissible. 
□ 

Step 2: construction of E tjk . Let £ = {a, S, ^ 2 } be a strongly 
admissible quadruple, and assume that i ^ a (a). Let £ t be the quadrupple 
defined by (4.21). Hence, according to the previous lemma, the quadruple £ 4 
is also strongly admissible for each tel. But then we can apply the results 
of Section 3 to associate with each £ t a family of quadruples E tjk , k G N 0 . 

Since E t is strongly admissible, we know from Lemma 3.2 that a is in¬ 
vertible. In fact, by item (ii) in Lemma A.l, all eigenvalues of a are in C+. 
Following the constructions given in Section 3, put 

#i(t, k) = ( I n + ia -1 ) fc -di(f), #2,k(t) = (4 - ia~ l ) k d 2 {t), (4.26) 

A(t, k) = [&i(t,k) d 2 {t,k)]. (4.27) 
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Furthermore, define S(t, k ) by 

S(t, k + 1) = S(t, k) + a^Sit, k)a~*+ 

+ a _1 S'(t, k)a~*A{t, k)jA(t, k)*, k = 0 , 1, 2 ,, 

S{t, 0) = S{t). 

Using these matrices we set 

E t) * = {a, S(t, k), (4 + foT 1 )^), (4 - icT 1 )*^)}. (4.28) 

By Lemma 3.2, the quadruple E tjk is also strongly admissible. 

Let E = {a, S, i?i, r? 2 } be a strongly admissible quadruple, and let 4 be 
the strongly admissible quadruple defined by (4.21). Following (3.12) put 

C k (t) =j+ A (t, k)*S{t, ky'Ait, k)+ 

— A(t, k + l)*S(t, k + l) _1 A(t, k + 1), k £ N 0 . 

Then, by Definition 3.3, the sequence { C k (t )} k€ ^ 0 is the pseudo-expontential 
potential generated by E f . The corresponding SkDDS is given by 

Vk+i(t, z) = (l m + z) (k G N 0 ). 

Furthermore, we know (see Propositon 3.6) that Ck(t) admits a factorization 
of the form: 


c k (t) = U k (t)jU k {t)\ U k (t)U k (ty = I m , ke No. (4.29) 

We also note that (4.9) and (4.10) are immediate consequences of the iden¬ 
tities (3.38) and (3.39). Indeed, since E t = {a, S{t), tfi(t), i4(£)} is strongly 
admissible, we just apply Lemma 3.10 with E t = E tj o in place of E 0 . Thus it 
remains to prove the identity (4.8) in order to prove Theorem 4.2. First, we 
formulate a proposition, which will be proved later. 

Proposition 4.5. Let E = {cc, So, 4, $ 2 } be a strongly admissible quadruple, 
assume thati ^ cr(a), and let E tjk be the strongly admissible quadruple defined 
by (4.28), where k E N 0 and tel. Then the function Y k (t , z) defined by 

Y k (t, z) = Ws tk (—z) ( 'l m + exp j— 2t (7^- + } > ( 4 - 30 ) 

where z y ±i and z a(—a), is a solution of equation (4.2) as well as of 
equation (4.3). 
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Proof of Theorem 4.2. As soon as Proposition 4.5 is proved we can prove 
Theorem 4.2 by using Proposition 4.1. Indeed, assume that the function Y 
defined by (4.30) satisfies (4.2) and (4.3). Then 

4-Yfc+i(t, z) = F k+1 (t , z)Y k+1 (t, z) = F k+ i(t, z)G k (t, z)Y k (t, z), 
at 

and 

^Y k+1 (t,z) = ^(G k (t,z)Y k (t,z) > j 

= (j t G k (t , z)^j Y k (t , z) + G k (t, z)^Y k (t, z) 

= (^G k (t,z)^ + G k (t,z)F k (t,z)'j Y k (t,z). 

We conclude that 

Gk(t , z )) Y k (t,z) = ^F k+1 (t,z)G k (t,z) -G k (t,z)F k (t,z)^Y k (t,z). 

Finally, note that z ^ Pi and z ^ a(—a) imply that Y k (t,z ) is non-singular. 
But then the above identity yields the compatibility equation (4.1), and we 
can apply Proposition 4.1 to get the identity (4.8). □ 

In order to prove Proposition 4.5 we need to compute the derivatives 

J f l\(t,k)’S(t,k)-\ jWv tk . (4.31) 

This will be done in a couple of steps. 

The first derivative in (4.31). Recall that A (t, k) is given by (4.26) and 
(4.27). It follows that 

(t, k) = [(I n + i(*)) (4 - ia ~ 1 ) k ■ 

But then, using (4.18) and (4.17), we see that 

= -2 [(a - 2 / n ) -1 (/ n + ict _1 ) fe i9i(i) ( a + il n ) -1 (/ n - ia~ 1 ) k d 2 (t)] 

= -2 [(a - il n )~ lr d i(f, k) (a + il n y l d 2 (t 1 k)] 

= -2 ((a - iQ-'Afa k)P\ + (a + k)P 2 ). (4.32) 
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Here P\ and P 2 are the projctions defined in (4.19) and (4.20), respectively. 

The second derivative in (4.31). In order to compute the second deriva¬ 
tive in (4.31) we present an alternative way to obtain the quadruple T, t ,k 
given by (4.28). As before, let E = {a, So, tfi, $ 2 } be a strongly admissible 
quadruple, and assume that i ^ a(a). We know that cr(ct) C C+. In a 
particular a is non-singular. As in Section 3, let 

E fc = {a, S k , ( I n + iaT 1 )*# 1 , (I n - ia~ l ) k d 2 }. 

From Lemma 3.2 we know that E*. is again a strongly admissible quadruple. 
The assumption i ^ a(a) allows us to apply Lemma 4.4 to E fc in place of E. 
Put 


tfi(M) = (I n + ia~ 1 ) k & 1 , (4.33) 

0 2 (M) = e~ 2t( - a+iI ^~\l n - ia (4.34) 
A (k,t)=[Mk,t) MKt)], (4.35) 

and let S(k,t ) be given by 

4s(m) = 

= — (a — i/ n ) -1 S(A;, t) + (a + iI n )~ 1 S{k, t) + 

+ S(k, t)(a* + il n )~ l + S(k, t)(a* - il n y l + (4.36) 

+ 2 (a 2 + Q- 1 (aA(k, t)jA{k, t)* + A (fc, t) jA(fc, t)*a*) ((a*) 2 + J n ) _1 ] , 
S{k,Q) = S k . 

Note that (4.36) is the analogue of (4.16) with E fc in place of E. Thus, by 
Lemma 4.4, the quadruple 

E k ,t = {a,S(k,t),tf i(/c,i),$ 2 (M)} 

is strongly admissible. 
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Lemma 4.6. For each k e No and t e R. we have T, t ^ = E^. In particular, 

j t S{t ' k) = 

= — (a — ■i/ ri ) _1 S'(t, fc) + (a + H n )~ 1 S(t, k)+ 

+ S(t, k)(a* + il n )~ 1 + A;) (a* - *7 n ) _1 + (4.37) 

+ 2(a 2 + 7 n ) _1 (aA(f, k)jA(t, k)* + A(f, k)jA(t, k)*a*^j ((a*) 2 + 7 n ) _1 . 

Proof. From (4.18) and (4.26) and the identities (4.33) and (4.34) we see 
that 

A(t, k ) = tfi(M) = e - 2t( “- iJ ") _1 (/ n + i£K- 1 ) fc i?i, 

tf 2 (*, *0 = tf 2 (M) = e- 2t ( Q+i7 ") _1 (/ n -za" 1 )*^. 

It follows that A(t,k) = A(k,t). Recall that E^ and E*.^ are both strongly 
admissible. Then item (ii) in Lemma A.l tells us that S(t,k and S(k,t ) are 
uniquely determined by A(t, k) and A (k,t). Thus S(t, k and S(k,t ) coincide, 
and hence E tjfe = E fc>t . Thus E t;fe = E fcjt . 

Finally, using A (t,k) = A(k,t) and S(t,k ) = S(k,t), we see that (4.37) 
follows from (4.36). □ 

The third derivative in (4.31). From (4.32) and (4.37) we shall derive 
the identity 

j t (A(t, k)*S(t, k)- 1 ) = 77+(t)A(t, k)*S(t , /c) _1 (a - i/ n ) _1 + 

+ 77 fe “(f)A(f, k)*S(t, k)~\a + il n ) _1 . (4.38) 

Here 77^ and 77^7 are the functions given by (4.19) and (4.20), respectively. 
To do this, note that according to (4.32) we have 

j t (A(t, k)*S(t, k)- 1 ) = -2PiA(f, A;)*(a* + i7 n )- 1 + 

-2P 2 A(t,fc)*(a*-27 n )- 1 + 

- A(t, A;)*S(t, k)- 1 ( j t S(t , A:))S(f, A:) -1 . (4.39) 

Identity (3.7) implies that 

( a* ± z7 n )^ 1 S'(t, A;) -1 = S'(7, fc) -1 (a ± i7 n ) _1 + 

+ i(a* ± f7 n )^ 1 5'(t, A;) _1 A(£, k)A(t, k)*S(t, k)~ x (a ± il n )~ l ■ (4.40) 
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By using (4.37), (4.40) and the equalities 

2(a 2 + InY 1 = i{(a + il n )~ l - (a -i/ n ) _1 ), 

2(a 2 + /n) -1 a = (a - i/ n ) _1 + (a + i/ n ) _1 , (4.41) 


after some calculations, we see that (4.39) can be rewritten in the form: 
j t (A(t,k)*S(t,k)-') = 

= V+(t)A{t, k)*S(t, k)~\a - tin)- 1 + V k ~(t)A(t, k)*S(t, k)~\a + iI n )-\ 

ffere 

v k + (£) = Im+j - iA(t, k)*S(t , k)~ 1 (a - iI n )~ l A(t , k)j+ 

+ ijA(t, k)*(a* — ■i/ n ) _1 S'(f, /c) _1 A(£, k) + A(£, k)*S(t , /c) _1 x 
x (a — iI n )- l A{t, k)jA(t, k)*(a* — iI n )- l S(t , /c) _1 A(f, /c), (4.42) 


V k (t) = I m -j+ iA(t,k)*S(t,k) 1 (a + il n ) 1 A(t, fc)j-h 

— ijA(t, k)*(a* + z/ n )^ 1 S'(f, /c) _1 A(t, k) — A (t, k)*S(t, /c) _1 x 

x (a + i/ n ) _1 A(f, k)jA(t, k)*(a* + i/ n ) _1 S'(f, fc) _1 A(f, /c). (4.43) 


It remains to show that I4 + (f) = H k {t) an d V). (£) = H k (t). To do this note 
that (3.51) and (3.52) tell us that H k (t) and H k {t) are also given by 


i/+(t) = I m + Wv t ' k (i)jWv t ' k (-iY, (4-44) 

tffc (*) = Im ~ W^ t k {-i)jW^ t k {iy. (4.45) 

Now use (4.47) and compute the products. This shows that H k {t) is equal 
to the right-hand side of (4.42) and H k (t ) is equal to the right-hand side of 
(4.43). In other words, V k + (t) = H k {t ) and V k (t) = H k (t), and hence (4.38) 
is proved. 

The fourth derivative in (4.31). In this part we shall show that 



UP 

A H - i 




2WV,(A) 


't,k 

Pi 

A — i 


: + 



(4.46) 
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(4.47) 


First, recall that W^ tk is given by 

W Etjfc (z) = I m + iA(t, k)*S(t, k)~ l (zl n - a) _1 A(t, k). 

Using (4.32) and (4.38), we obtain 

(^Ws t , fc ) (\)=i(H+(t)A(t,kyS(t,k)- 1 (a-iI n )- 1 {\I n -a)~ 1 A(t,k)+ 

+ H^(t)A(t, k)*S(t, ky l (a + i/ n ) _1 (A/ n - o) _1 A(f, /c)+ 

- 2A(t, fc)*S(t, - il n )-\\l n - a)-^^, fc)Pi+ 

- 2A(t, fc)*S(t, fc) _1 (a + i/ n ) _1 (A/ n - a)"^, fc)P 2 ). (4.48) 

Clearly we have 

(a - H n )~ 1 (XI n - a) -1 = (A - i) _1 ((a - f/„) _1 + (A/ n - a) -1 ), (4.49) 

(a + i/n) -1 (A/„. - a) -1 = (A + i) _1 ((a: + U n )~ l + (A/ n - a) -1 )- (4.50) 

Taking into account (4.47), (4.49) and (4.50), we rewrite (4.48) as 

W = (A - (wy.W - ny,M)+ 

+ (A + i)- 1 ff,7(i)(ir E .,(A) - IU, ,,(->)) + 

- 2(A - y A)-Ws, ,(;))/>, 

-2(A + i)- 1 (Vr ElJ ,(A) - WVy-0)A. (4.51) 

Applying (3.46) and (3.47) with E^ in place of E/, ; we see that 

= 2 W Eu (i)P u H^(t)Wz lt (-i) = 2W Ell (-i)P 2 . (4.52) 

In view of (4.52), formula (4.51) can be rewritten as (4.46). 

Proof of Proposition 4.5. Let Y be the matrix fucntion defined by (4.30). 
The fact that Yk(t,z ) is a solution of equation (4.2) follows by applying 
Theorem 3.5 with S 4 in place of E. Similarly, using (4.46) and the fact that 
the matrices Pi, P 2 , and j are mutually commutative, we see that Y k (t, z) is 
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a solution of the system (4.3). Indeed, 


d 


dt 


d 


—Y k (t, z ) = — We, fc (--) Im + -j exp -2 1 —— + 


Pi_ P 2 


dt"^\ 

-2 Y k (t,z) ^ z+ . , . 

+ Yk (t, z) - 2 Y k (t, z) 

—z — l —Z + lJ 

p i_ + P 2 


P 1 , P 2 


Z + l z — i 


p i _ p 2 


+ 


—z — i —z + i 


- 2 Y k (t,z) 


z+i z—i 


= _(ym + sm YtM . 

\ z +« z — l J 
Thus Y has the desired properties. 


□ 


Remark 4.7. The equality H k = 2I m — (H k )* is immediate from (4.44) and 
(4.45). 


A remark on the continuous space variable x analogue. The transfer 
matrix function of the form (2.10) depending on the continuous parameter 
t and discrete parameter k, which we use in this section, appeared in [27], 
whereas the transfer matrix function (of the form (2.10)) depending on con¬ 
tinuous parameters x and t appeared first in [36,37] (see also [22] and [23]). 
To explain the case of two continuous parameters in more detail we begin 
with the following analogue of Lemma 4.4. 


Lemma 4.8. Let £ = {a, So,0i,0 2 } be a strongly admissible quadruple, and 
let k be a positive integer. Define 


Mt) = e~ itak Vi, 0 2 (t) = Ht) = [0i(t) Mt)] , 

t k 

S(t) = Sq + f '^^a u ~ 1 A(s)jA(s)(a*) k ~ u ds, t > 0. (4.53) 


0 V=1 


Then for each 1 6 1 the quadruple £(i) = {a, S(t), flfit), 0 2 (^)} Is strongly 
admissible. 


The proof of the above lemma follows the same line of reasoning as used 
in the proof of Lemma 4.4. We omit the details. 
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Next, using the above lemma, we construct a family of pseudo-exponential 
potentials depending continuously on an additional variable 1 6 i. The 
starting point is the pseudo-exponential potential v in (2.1) defined by the 
strongly admissible quadruple £ = {a, S'o,'di, ^ 2 }- Fix i G R, and let £(t) 
be the strongly admissible quadruple defined in Lemma 4.8, i.e., £(t) = 
{a, S(t), 'di(t), $ 2 (£)}- In particular, S(t) is given by (4.53). We apply 
Lemma 2.1 with £(f) in place of £. Put 

A(x,t) = e~ ixa A(t) = 

tfaOM) = e ixa ®2 (t) = e i{ - xa+itak) d 2) 

A (x,t) = [fii(x,t) $ 2 (M)] , 

S(x,t) = S(t) + / A(s,t)jA(s,t)*ds, x>0. (4.54) 

Jo 

Then, by Lemma 2.1, the quadruple 

£(x,f) = {a, S(x, t), $i(x, t), t)} (4.55) 

is admissible for each x > 0. Furthermore, the corresponding pseudo¬ 
exponential is given by 

v(x, t) = 2$* 1 e i ( xa * +it ( a * )k) S(x, x > 0 (4.56) 


In this way we obtain a family of pseudo-exponential potentials depending 
on the additional parameter t. 

For k = 2,3 the potential v in (4.56) satisfies important integrable nonlin¬ 
ear equations; see, e.g., [36, Section 2] and [23, Section 4], Namely, for k = 2 
the function v is a solution of the matrix nonlinear Schrodinger equation 


dv d 2 v * 

2 di + t aj + 2ivvv = (> ’ 


(4.57) 


and for k = 3 it satisfies the matrix modified Korteweg-de Vries equation 


dv d 3 v 

4-)- 

dt dx 3 


+ 3 


dv 


;dv\ 


„ + „ , 

ox dx/ 


= 0 . 


(4.58) 


The proofs of the above results can be obtained by direct computa¬ 
tions as in [36, Section 2] and [23, Section 4] or by applying the generalized 
Backlund-Darboux transformation [38] to auxiliary linear systems and using 
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compatibility condition (zero curvature equation) G t — F x + GF — FG = 0, 
which is a continuous case equivalent of the compatibility condition (4.1). 
For instance, equation (4.57) is equivalent to the compatibility condition 
Gt — F x + GF — FG = 0, where G = izj + V ( x , t ), 

F = i(z 2 j — izjV(x,t) - (V x (x,t)+jV(x,t) 2 )/ 2), 
and V has the form (1.2). 


A Appendix: admissible quadruples 

Let Sq and a be n x n matrices, and let and $2 be matrices of sizes n x mi 
and n x m 2 , respectively. The quadruple {a, Sq, $ 1 , $ 2 } is called admissible 


(see the beginning of Section 2) if 

Sq > 0 and aS 0 — S 0 a* = i AA*, where A := [$1 $ 2 ]- (A.l) 

By definition the transfer function associated with the quadruple 

£ = {a, S 0 , $2} 

is the matrix function We given by (see (2.10) and (3.15)): 

We(z) — I m + *A *Sq 1 (zI ti — a;) _1 A, where m := m 1 + m 2 . (A.2) 

With £ we associate two other rational matrix functions: 

ViM*) = i’d* 2 SQ 1 {zI n - fdi)~ 1 'd 1 , where := a - (A.3) 

if 2 ,x(z) = -ifiis^izln + where /3 2 := a - id^S q" 1 . (A.4) 

If £ = {a, S'o,'di, ^ 2 } is an admissible quadruple, then the same holds true 
for the associate quadruple 

£ # :={a,S o ,0 2> M (A.5) 

and a direct computation shows that 

7h,s#(--) = ¥>2 ,e(z)- (A.6) 
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Admissible quadruples are closely related to symmetric S-nodes (see [44, 
Section 2.2]). Indeed, if S = {a, S 0 , $ 1 , $ 2 } is an admissible quadruple, then 
the matrices 


A := a, S := S 0 , IU : = A, I\ := Sq% 

form a symmetric S-node (with J = I n ) as defined in Section 2.2 of [44], Using 
this connection, it is readily seen that item (iii) in the following lemma is a 
special case of Theorem 2.2.1 in [44] (see also [42], Theorem 17.1 in [7] or 
Corollary 1.15 in [41]). 

Lemma A.l. Let E = {a, So, $ 1 , $ 2 } be an admissible quadruple. Then 

(i) a (a) C C+; 

(ii) if, in addition, the pair {a, A} is controllable, then cr(a) C C + , and in 
that case S 0 is uniquely determined by a and A. More precisely, S 0 is 
given by 

1 f°° 

So = — / (A I n - a) _1 AA*(A/ n - a*y l d\. (A.7) 
2vr J„ 00 

(iii) W E is a rational Schur class function and its values on the real line are 
unitary matrices. Moreover, we have 

Im - W E (z)W E (C)* = 

= i(z - C)A *So 1 (zI n - a)- 1 So(C4 - a)- 1 S 0 ” 1 A, (A.8) 

Im - W E (z)*W E (z) = i(z - z)A*(zI n - a*)~ 1 sy(zl n - a) _1 A. (A.9) 

Proof. In remains to prove items (i) and (ii). Since So is positive definite, 
the second part of (A.l) implies that 

(S" 1/2 aS 0 1/2 ) - (S 0 -1/2 aS” 1/2 )* = zS" 1/2 AA*S” 1/2 , S" 1/2 AA*S" 1/2 > 0. 

It follows that the numerical range of the matrix Sq'^oSq' 72 is a subset of 
C + , and hence cr(a) = cr(S 0 " 1 ^ 2 aSo^ 2 ) C C+. This proves item (i). 

The hrst statement in item (ii) is a straightforward application of the 
classical Chen-Wimmer inertia theorem (see Corollary 1 in [33, Section 13.1]). 
For (A.7) see [26] or Theorem 1.4.1 in [21] (see, also (3.10)). 
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Lemma A. 2. Let £ = {a, So, i?i, $ 2 } be an admissible quadruple, and parti¬ 
tion its transfer function Wy as a 2x2 block matrix function, 


Wyfz) = 


, where a(z ) has size m\ x m\. 


(A.10) 


(A.ll) 


a(z) b(z) 
c(z ) d(z) 

Then the associated functions (pip and tp2 ,s are given by 

(pi ,x(z) = c(z)a(z )~ 1 and (p2,x(z) = &(— 2 )d(—z ) -1 

Proof. From (A.2) and A = [i?i $ 2 ] it follows that 

a(z) = I mi + id\Sf l (zI n - 0(2:) = 

Since a is a proper rational matrix function whose value at infinity is the 
m-] x m-] identity matrix, a classical inversion theorem (see, e.g., Theorem 
2.2 in [6]) tells us that 

a(z )" 1 = I mi — i'&fSf 1 (zl n — cx x ) _1 'di, where a x = a — ididlSf 1 = fd\. 

Hence, again using Theorem 2.2 in [6] and taking into account (A.3), we 
know that 


c(z)a(z) 1 = id* 2 S 0 l (zI n — a) 1 , d 1 a(z) 1 = id* 2 S Q 1 (zl n — a x ) 1 'd i 
= id^S^izIn - ff = ipi t y(z), 


which proves the first part of (A.ll). 

To prove the second identity in (A.ll) we use the associate admissible 
quadruple £ # = {a, S 0 , $ 2 , $ 1 } and equality (A.6). Note that 


\p2 $ 1 ] = [$1 $ 2 ] U, where U = 

Since U is unitary, it is clear from (A. 12) and (A.2) that 1F S # 
But then, using the special form of U, we obtain 


0 I mi 
Im 2 0 


(A.12) 


U*WyU. 


Wy# = U*WyU 


d(z) c(z ) 
b(z ) a(z) 


(A.13) 


Applying the result of the first part of the proof to £# in place of £ we see 
that <£>i )S # = b(z)d(z)~ 1 . But then (A.6) yields the second identity in (A.ll). 
□ 
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Theorem A.3. Let p be an m 2 x m\ rational matrix function. Then p 
coincides with some function p\p associated with an admissible quadruple £ 
if and only if p is strictly proper. In that case, the corresponding £ can be 
obtained explicitly by using the following procedure. 

Step 1. Let n be the McMillan degree of p, and construct a minimal real¬ 
ization of p : 

p(z) = id 2 (zl n - 7) _1 0i. (A. 14) 

Step 2. Choose X to be the unique positive definite solution of the algebraic 
Riccati equation 


7 x - x 7 * - ixe 2 e* 2 x + io^i = o. (A. 15 ) 

Step 3. Put 

So = In, $1 = x~ 1 / 2 0 u d 2 = x 1 ' 2 d 2 , (A. 16 ) 

a = A' 1/2 7 A 1/2 + mMJ. (A.17) 

With these choices the quadruple £ = {a, S 0 , $ 1 , $ 2 } is admissible and p 
coincides with its first associate function Moreover, in this case the 

pairs {a, D 1 } and {a,d 2 } are controllable. 

Proof. First note that the minimality of the realization (A.14) is equivalent 
to the requirement that simultaneously the pair { 6 2 , 7 } is observable and the 
pair { 7 , 6 1 !} is controllable. But then we can apply Proposition 2.2 in [23] to 
show that the Riccati equation (A.15) has a positive definite solution X, a 
result which has its roots in Kalman’s theory of mathematical systems [28] 
(see also [32, pp. 358 and 369], where the uniqueness of X is shown as well). 
The remaining part of the proof is split into two parts. 

Part 1. Using the definition of a in (A. 17), the unique positive definite 
solution X of the Riccati equation (A.15), and the definitions of 'd\ and d 2 
in (A. 16) we see that 

a - a* = A“ 1/2 7 A 1/2 - X 1/2 7 *X" 1/2 + 

= ix 1 / 2 e 2 e* 2 x 1/2 - •la” 1 / 2 m*x - 1/2 + 2^1 

= iD 2 $2 — i$i$l + 2iti 1 -&\ = HI 2^2 + i'&i'dl- 
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Since S 0 = I n by definition, we conclude that £ = {a, S 0 , i9i, d 2 } is an 
admissible quadruple. Moreover, the associate function <p 1 is equal to p. 
Indeed, in this case using (A. 17), we have 

Pi = a- itiiVl = X~ 1/2 ^X 1/2 . (A.18) 

Thus X 1 / 2 /3iX -1//2 = 7, and hence 

<Pi,z {z) = i’d* 2 SQ 1 {zI n - Pi)~ l di 

= ie* 2 x 1 / 2 (zi n - p l )- 1 x- 1 / 2 e 1 
= i0 2 (zi n - x l/ 2 p 1 x- l/ 2 )~ 1 e 1 

= i0 2 (zl n - 7 ) _1 6'i = ip. 


Part 2. In this part we prove that the pairs {a, $ 1 } and {a, $2} are control¬ 
lable. Since the realization (A. 14) is minimal and p = pip, the same is valid 
for the realization (A.3), and hence the pairs {/Si, $ 1 } and {(/fi)*, $ 2 } are con¬ 
trollable. Formula (A. 3 ) implies that a — P\ -\-id\d\. Using the latter iden¬ 
tity, the fact that {/Si,i9i} is controllable and the equality lin'd! = Im^i?*, 
we conclude that {a, $ 1 } is controllable too. 

It remains to prove that {a,i9 2 } is controllable. Since S 0 = I n and the 
quadruple £ = {a, Sq, i9i, i 9 2 } is admissible, we have a — a* = i(didl + d 2 dlf), 
and therefore 


Pi = a — id = a* + id 2 d 2 = (a — id 2 d 2 )*. 

Thus a = (Pi)* + id 2 d 2 . We know that {(/Si)*, i 9 2 } is controllable and that 
Imi9 2 = 1111192192. These facts imply that {a,i9 2 } is also controllable. □ 

Using the identity (A.6) and Theorem A.3 we can show that a proper 
rational matrix function p also coincides with the second associate function 
<^ 2 ,s of some admissible quadruple £. More precisely, we have the following 
corollary. 

Corollary A.4. Let p be an m 2 x mi rational matrix function. Then p 
coincides with some function p 2 pi associated with an admissible quadruple £ 
if and only if p is strictly proper. The corresponding £ = {a, S' 0 ,i9i,i9 2 } can 
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be obtained explicitly by using the following procedure. First, we construct a 
minimal realization of dp, 


y(z) =-i6l(zl n + 7 ) 1 9 2 . 


Next, we choose X to be the unique positive definite solution of the Riccati 
equation 'jX — X 7 * — iX9\9\X + i9 2 9\ = 0. Finally, we put 


So = In, 9i =x 1 ' 2 9 1 , ■& 2 =X- 1 ' 2 9 2 , 

a = X~ 1/2 ^X 1/2 + id 2 9* 2 - 


(A.19) 
(A.20) 


Then p = p 2 ^ and the pairs {a,$i} and {a,d 2 } are controllable. 

Proof. Put p{z) = p(—z). Then p is rational, and p is strictly proper if 
and only if p is strictly proper. But, by Theorem A.3, the latter happens 
if and only if there exists an admissible quadruple £ = {a, So, $ 1 , $ 2 } such 
that p(z) = p l j ] (z). Next, taking into account (A.5), we see that £ is the 

associate quadruple £ # for £ = {a, So, 9 2 , $ 1 }. Then we know from (A.6) 
that 

p( z ) = p(-z) = Pi )12 #{-z) = ^ 2 ,e0). 

Moreover, Theorem A.3 provides a method to construct £, which in terms of 
£ = £ # (or, equivalently, £ = £#) yields the procedure to recover £ from 
p and formula (A.4). □ 

Proof of Proposition 2.8. Let S(x) be the n x n matrix defined by 
(2.2), and put 


— 2ixA d n 

0 ’ 


Q(x) = [S 0 -U n \ e 




where A is given by (2.23). In order to prove Proposition 2.8, that is, to prove 
(2.24), it suffices to show that Q(x) = X(x), where A(x) = e~ lxa S(x)e~ lxa *, 
for each x G M + . For that we directly differentiate A(x) using the first 
identity in (2.7) (where we substitute the expression for A from (2.2)) and 
obtain 
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Now, we show that Q satisfies the same first order linear differential equation 


-^Q( x ) — —i(aQ(x) + Q(x)a*) + e 2lxa did\ - d 2 d* 2 e 
To prove (A.22) note that 

i-^Q(x ) = [S 0 -il n ] (Ae" 2 ^ + e" 2 ^) 

Using (2.3) and the definition of A in (2.23) we see that 


[So —H n \ A = [So«* + id\d\ —iot ] 

= [aS 0 — id 2^2 — ioi\ 

= a [S 0 -il n \ - i [d 2 d * 2 0 ] . 

Furthermore, we have 

•4 fo] = [_£*.] = [o] “* + [/„] 

It follows (using again the definiton of A in (2.23)) that 


—2 ixa* 


(A.22) 


i-^Q(x)=a[S 0 —il n ] e~ 2ixA ^ -i [d 2 d* 0] e" 2 “ A 


a 


+ [So —iln] e~ 2ixA 
+ [So -iI n ]e~ 2ixA " (-W1) 

J-n 

= aQ(x) - id 2 d* 2 e~ 2ixa * + Q(x)a* + ie-^d^. 


(A.23) 


Equation (A.22) is immediate from Q. Since A and Q satisfy the same 
first order linear differential equation with the same initial condition A(0) = 
Q(0) = S 0 , we see that A(x) = Q(x) for each x E R+, as desired. 
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